Ginzburg-Landau theory of dirty two band s± superconductors 
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In this paper we study the effect of non-magnetic impurities on two-band s± superconductors 
by deriving the corresponding Ginzburg-Landau (GL) equation. Depending on the strength of 
(impurity-induced) inter-band scattering we find that there are two distinctive regions where the 
superconductors behave very differently. In the strong impurity induced inter-band scattering regime 
Tc « where rt ~ mean-life time an electron stays in one band the two-band superconductor 
behaves as an effective one-band dirty superconductor. In the other limit Tc > Tf^, the dirty 
two-band superconductor is described by a network of frustrated two-band superconductor grains 
connected by Josepshon tunnelling junctions. We argue that most pnictide superconductors are in 
the later regime. 



With the discovery of the Iron-based (pnictides) su- 
perconductors, superconductivity characterized by more 
than one order parameters, i.e. the multi-gap supercon- 
ductors, becomes a topic of interests. Band structure 
calculations indicate that the materials have a quasi-two- 
dimensional electronic structure, with five bands centered 
around the T- and M- points in the Brillouin zone con- 
tributing to the Fermi surface. It has been proposed that 
the superconducting order parameters in this multi-band 
materials has so called s±-wave symmetry, where the or- 
der parameters have s-symmetry but with opposite sign 
between bands centered at T- and Af-points[l|, [3, 

The effect of impurities in this class of materials has 
been an issue of interests. NMR[4] and lower critical 
field dataQ seems to suggest the existence of nodes 
in the superconducting order parameter while APRES 
experiemntsP, ^ favor node-less gaps. One possible 
solution to this controversy is that large number of in 
gap states are induced by impurities in the material be- 
cause of the special s± order-parameters. Indeed, such a 
scenario has received supports from self-consistent-Born 
type calculation where in gap states are found to appear 
easily in s± superconductors [S, t^. 

In this paper we study the effect of non-magnetic im- 
purities on two-band s±-wave superconductors by ana- 
lyzing the corresponding Ginzburg-Landau theory. The 
effect of impurities is included by generalizing the stan- 
dard Bogoliubov de Gennes theory and diagram- 
matic perturbation techniques 11| which have been ap- 
plied to study the effect of impurities on single-band 
superconductors 13, 13, l3l to the case of two-band su- 
perconductors. 

We start with the Bogoliubov-de Gennes formulation 
of BCS theory ff^. The system we consider is charac- 
terized by a BCS Hamiltonian, H = Hq + Vbcs, where 



i,j,(T 

where i,j = 1,2 and <J are the band and spin in- 

dices, respectively. Hoi{V) is the band Hamiltonian de- 
scribing electronic wave-functions in band i. Uij{r) is 



a non-magnetic disordered potential which scatters elec- 
trons both within [i — j) and between (z ^ j) bands. 
ipjai'4'j'ij) electron annihilation (creation) operators. 



Vbcs = - E / d''r^tA^i^tAr)i'M^i'jAr), 



(lb) 

is the BCS interaction between electrons, a — —a. We 
note that Vij > means attractive interaction in our 
notation. Introducing the BCS decoupling. 



where = {ipid-4'iu) we obtain the Bogoliubov-de 
Gennes equations for quasi-particle states n[lO|, 

enul:H^ = 5](%i?o.(V) + C/.,(f))zi(^)+A,(r>«(f;2a) 



en^^W (^1 = - E ('^'^■^0*. (V) + U., (r)) vlP + A* (f),.W (r) 



where Ai(r)'s are determined by the self-consistent equa- 
tion, 

A,(f) = - ^ V,,A,(r} = V.,ui^\f)vi^>{f){l - 2/„), 

j j,n 

(2b) 

where /„ = \/{e^^^'^" + 1) is the Fermi-Dirac distribu- 
tion. 

We note that inter-band electron pairing is not in- 
cluded in our mean-field BCS decoupling. Physically 
different electronic bands describe electrons located at 
different parts of the Brillouin zone and an inter-band 
pairing implies finite-momentum Cooper pairs which is 
usually energetically not favorable. The mean-field de- 
coupling we employed introduces only Josepshon cou- 
pling between superconducting order parameters in the 
two bands and the electronic wave- functions in the two 
bands are mixed only by the disorder-potential Uij. 

The Ginzburg-Landau (GL) equation for the system 
can be derived by assuming that Ai(r) is small and ex- 
panding Eq. ([2]) in powers of Ai(r) to third order. We 
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furthermore assume that Ai{r) is slowly varying and 
perform a gradient expansion Ai(f') ~ Ai{f) + (r' — 
r).VAi(f) + ... to obtain 

A,(f) = ^fe°)(rOA,(r) + iifW(rOV2A,-(f))(3) 



-5]4"LWA,(f)A^(f)A,(f), 



dV(f-r')2"if,,(r,r'), (4) 
d'^n / d'^rz / d'^r3L,jki{r,fi,r2,r3), 



where 
where 

k iujn 

(5a) 

Lijki{ri,r2,f3,rA) = - ^ V^^m x ^ ^ 5,„j(fi, f2, <<&b) 
gmk{ri,rz, -iuJn)gi]ir4,r2,iuJn)gikiri,f3, -iUn) 
and 



5y 



(5c) 



(l>n\f)'s are eigenstates of i/o given by en(f)n' (r 

Ej (%i?o»(V) + [/.,(f)) #(f). 

To study the effect of impurities we first consider 
the impurity-averaged GL equation where we replace 
K<-°\K'^^^ andL(o) by their averages over disorder poten- 
tial Uij{f). Notice that we have assumed that {KA)av ^ 
{K)avA, etc. in this process where {.■)av denotes impu- 
rity average [l^ . The validity of this approximation 
will be examined later. With this approximation we ob- 
tain the usual impurity-averaged GL equation 

A.(r) = (KS^^jir) + ii^i;V2A,(f)) (6) 



jki 

where X*^") = (/^("^(r))^!,, etc. We shall consider the 
limit EfT >> 1 where Ef ^ Fermi energy and r ^ elastic 
scattering life time in our calculation and compute the 
impurity average to lowest order in impurity density rii 
(semi-classical limit) [13, 12, 13|. In this limit electron 
motion becomes diffusive and modifies the long-distance 
behavior of if ^"^ and L^^li . 



To compute {Kij{f,f'))av we note that it can be writ- 
ten aslisl 



where 
T { 

^m)) av 



'))o 
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dEdE' 



function of 
the Fourier 



Fk,{E - E',r 
4'^(r)<^i^>(r')4^Hr)0(^)' 
is related to the 
the 



FkjiE-E',r-f) 
{iujn - E){-iu>n - E') 
(7) 



{r')5{E - e^)5{E' - 
to the density-density response 
corresponding dirty metal [lH. 



= Imxij{q,u} + iS), where Fij{uj,q) is 
transform of Fij[uj,r) and Xiji^i^ + *'^) 
is the (zj) component of the density-density response 



function of the dirty metal [111. Ll 



The density-density response function can be evalu- 
ated to lowest order in impurity concentration by keep- 
ing the lowest order self-energy and particle-hole ladder 
diagrams [in. To perform the impurity average we as- 
sume {Ui-i{r))av = and {Ui.j{r)Uki{r'))av ^ only if 
i = k,j = I or i = l,j = k with {Uu{r}Uti{f'))av = 
<5'*(f-r')n,|u,|2 and {U,2(2i){f^U,2(2i)ir'))av = S^if - 
r')ni\ut\'^, i.e. the different type of scattering events are 
uncorrelated with each other. The corresponding aver- 
aged retarded (R) and advanced (A) electron Green's 
functions have the form 11 1 



{g^^^^\k,Uj))av 



-^^k + i-)^ 



where ^ = r„; ^ + t^j\ r,/ = 2niTr\ui\^ N,{0) and 
T^^ = 2nin\ut\'^N-i{0), where 1(2) = 2(1) and Ni{0) is 
the density of states for band i electrons on the Fermi 
surface, is the mean life time where an electron in a 
state in band i is scattered to another state in band j. 
Notice that the impurity-averaged Green's function has 
no off-diagonal (z ^ j) term. 

The corresponding density-density response function 
is calculated to lowest order in rii by summing ladder 
diagrams in particle-hole channel (fig.l). We shall be 
interested at the low energy, long wave-length transport 
behaviors of the system. In this limit we need to keep 
only those processes where the particles and holes are 
coming from the same band in our calculation. This is 
because the two bands are located at different parts of 
the Brillouin zone, and the center of mass momentum of 
inter-band particle-hole excitations are usually large and 
do not contribute to small q processes. 

Evaluating the diagrams, we obtain 



Xii{q,^) 



and 



{-iuj + D2q^ + T2i)i-iuJ + + rf/) - (riarsi)"! 

(8a) 



Xiiiq,^) = 



{-ioj + D2q^ + T2i){-ioJ + Diq^ + r^/) - (narsi)-! 



(8b) 
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Figure 1: ladder diagrams in particle-hole channel, i,j,k are 
band indices. We include only processes where particle and 
hole are coming from the same band. 



where _D,j — VpTi/d is the diffusion constant for band 
electrons and 



=P,;(0)P,(0)7V,(0), Pi(2)(0) 



^1(2) (0) 

Nt{0) 



(9) 



where 7Vt(0) = 7Vi(0) + A^2(0) is the total density of states 
on the Fermi surface. The result is valid in the small q, lo 
limit uj « and Ag^ << for both i = 1, 2. 

Fij{q,Lu) can be evaluated using Eq. ([8]) and has very 
different behaviors at energy scales higher and lower than 
the inter-band scattering life-time t~^. To simplify cal- 
culation we shall assume T12 ~ T21 ~ are of the same 
order of magnitude. In this case we obtain in the limit 
UJ « T^r^ and Diq"^ « rj"\ 



Frj{q,uj) 



P^jDeffq^ 
LO-^ + {D,ffq^)^' 



(10a) 



where A// = Pi(0)i:'i + P2(0)P'2 and 



F,j{q,uj)^S,jN,{0)- 



D,q' 



{D^q^r'' 



(10b) 



in the opposite limit uj >> t^^ and Diq^ >> rf^. Phys- 
ically, electrons have scattered many times between the 
two band already in the limit iu,Diq^ « rf^ and the 
identity of bands is lost as far as electron dynamics is 
concerned. The only remaining information of "bands" is 
that electrons have probability Pj(0) of residing in band i. 
The identity of the two bands remain in the opposite limit 
UJ, Diq^ >> T^^ where electrons stay mainly in one-band. 
The two different limits expressed themselves in the GL 
equation where we find that in the limit Tc << t^"^, elec- 
trons have to scatter between the two bands many times 
before forming a Cooper pair and the identity of intra- 
band Cooper pairs is lost, whereas intra-band Cooper 
pairs survived in the opposite limit r^"^ << Tc- We shall 
first consider the limit Tc « rf^ in the following. 



Putting together Eqs. (U), (O and ([TO]) , we obtain 



[ID 



K. 



(1) 



Vik X pkj — — 

k 



where ujd is the high energy cutoff for the attractive in- 
teraction in BCS theory. We have assumed t^^ » Tc 
in deriving . 

-^ijki '^^^ computed similarly in perturbation theory. 
We obtain after some lengthy algebra 



.(0) 



^ ^ V„nPrn{0)Pj{0)Pk{0)Pl{0)Nt{0). (12) 



The result can be understood most easily by noting that 
in the limit Tc << rjT^ the dynamics of electron is de- 
scribed in an effective single-band picture with probabil- 
ity Pi(0) of finding electrons in band i. 

Putting Eqs. (g]), ^ and ^ in Eq. (O, multiple 
the resulting equation by Pi(0) and sum over z, we obtain 
an effective single band GL equation 



Acffir) = a{T)Acff{r)^-V'Acff{r^-c\Acff{r)\'Acffir) 

(13) 

where A,//(7=) = P,(0)A,(7=) and a(T) - 

K.A^t(0)ln(wd/P), b ^ VaMQ)Dcff/Tc, c ^ 
VavNtiO)/T^ where = E^J P^iO)V,jPJ{0) is the aver- 
age interaction electrons see in forming the Cooper pairs. 

The individual band order parameters Aiif) are re- 
lated to Acff{r) by 



A.(r-) 



1 



{Y,VMO))Acff{f). 



(14) 



and are 'slaved' to Ae// in the sense that they are not in- 
dependent dynamical variables in the system. The dirty 
two-band superconductor behaves as an effective dirty 
one-band superconductor in the regime Tc << t^^ where 
measurement of superfluid properties cannot distinguish 
between whether the system is originally a single-band 
or a two-band superconductor. 

The effective single-band description has a number of 
interesting predictions. The (average) superconducting 
transition temperature is given by 



: uJd exp( 



(15) 



which is very different from clean two-band superconduc- 
tors where Tc is determined by 



(16) 
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where Vu 

Vt = Vi2 

sgn{Vt). 



VuN.iO) and Vt = Vt^ Ni{0)N2{0) where 
V21. Notice that ri"'' is independent of 



It is straightforward to show that Tc < Tc^\ i.e. Tc is 
always lowered by disorder. However Eq. pS]) says that 
the precise value of Tc is insensitive to the strength of dis- 
order and depends only on the density of states of the two 
Fermi surfaces in the limit t^T^ << TJ This surpri sing 
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result is a direct consequence of "Anderson Theorem 
applied to the (effective) one-band superconductor. 

Contrary to the case of clean superconductors we also 
observe that Tc depends now on the sign of Vt . In partic- 
ular Tc is enhanced by Vt only if Vt > 0, suggesting that 
disorder disfavor s± state. The relative sign between Ai 
and A2 depends on all the interactions now (Eq. (fT4|) ) 
and is not solely determined by sgn(Vt)! 

Next we consider the regime t^^ < Tc- This region is 
non-trivial as can be seen from the change in Tc as a func- 
tion of Tj~^ determined by the GL theory. At r^^^ —^0 



Tc is determined by Eq. p6)) for clean superconductors 
whereas Tc is determined by Eq. (fTS]) at Tj~^ » Tc- Tc is 
different but insensitive to disorder at both regimes (An- 
derson Theorem)! Therefore Anderson Theorem must 
breaks down and Tc becomes sensitive to disorder at the 
intermediate regime < < Tc. The non-trivial effect 
of impurity scattering in this regime is shown in single- 
impurity calculations where it is found that in-gap bound 
states are induced easily by inter-band impurity scatter- 
ing and the Josephson coupling between the bands is sup- 
pressed correspondingly in the s± stated, S EBl- 
note that the in-gap states are absent in the rjT^ >> Tc 
limit where an effective single-band description becomes 
valid, consistent with findings on superconductors with 
sign-changing order-parameters [lij . 

The rare (but strong) effects of inter-band impurity 
scattering suggests that the self-averaging approximation 
{KA)av ~ {K)avA breaks down in the regime r^"^ < Tc 
and Ai{r) becomes sensitive to the precise configuration 
of inter-band scattering potentials. The sensitivity of 
Ai(r) to the impurity potential can also be seen directly 
from the (averaged) GL equation. It is easy to show that 



K, 



(1) 



-V,N,{Qi)D,/Tc^ 



and the GL equation does not take the form of an ef- 
fective single-band GL equation in this regime. As a 
result its solutions are very sensitive to local variations 
ini^f (r). 

Therefore to describe the effects of order at this regime 
we should start with the un-averaged equation ((3|). It is 
more convenient is to replace the continuum GL equation 
by a random Josephson coupling lattice model with free 



energy 

F = ^ [^a„,z(r;*)(A+(*)AK*)+c.c.) + ^6„,|A,„(z)(l^) 

I \ m/ rn / 

<ij>,i 

where {i,j) and {l,rn) are lattice site and band indices, 
respectively. < i,j > denotes nearest neighbor pair sites. 
The first term in (fT7|l represents grains of two-band su- 
perconductors where the two bands are coupled only 
through Josephson coupling 012(21)- The second term 
represents Josephson coupling between nearest neighbor 
grains. a;m(T;i) aim{T) with 012(21) > for clean 
s± superconductors and aim becomes randomized in the 
presence of disorder. It is easy to see from a three-site 
calculation that the phase of the order parameters are 
frustrated if sgn{ai2(2i){T]i)) becomes randomized fl7| . 
indicating that a uniform superconducting state becomes 
unstable when inter-band impurity scattering is strong 
enough. 

Experimentally, we note that different superconduct- 
ing gaps were observed at energy bands located at the 
T— and M— points of the pnictide superconductors 
in ARPES experiments Q, indicating that the materi- 
als are located in the weak inter-band scattering regime 
T(~^ < Tc where impurity- induced in-gap bound states 
are present, consistent with the existence of large den- 
sity of in-gap states found in NMR0] and lower critical 
field Q experiments. We propose here that a uniform su- 
perconducting state may become unstable at this regime. 
A detailed analysis of the superconducting behavior at 
this regime will be the subject of a separate paper. 
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